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Abstract 

We introduce a simple generalization of rational bubble models which removes the funda- 
mental problem discovered by [[20| that the distribution of returns is a power law with exponent 
less than 1, in contradiction with empirical data. The idea is that the price fluctuations asso- 
ciated with bubbles must on average grow with the mean market return r. When r is larger 
than the discount rate r$, the distribution of returns of the observable price, sum of the bubble 
component and of the fundamental price, exhibits an intermediate tail with an exponent which 
can be larger than 1. This regime r > r$ corresponds to a generalization of the rational bub- 
ble model in which the fundamental price is no more given by the discounted value of future 
dividends. We explain how this is possible. Our model predicts that, the higher is the mar- 
ket remuneration r above the discount rate, the larger is the power law exponent and thus the 
thinner is the tail of the distribution of price returns. 



1 The fundamental constraint on distribution of returns of rational 
bubbles 



Since the publication of the original contributions on rational expectations (RE) bubbles by [^] 
and [fflj, a huge literature has emerged on theoretical refinements of the original concept and the 
empirical detectability of RE bubbles in financial data (see [f7|] and [jl|], for surveys of this literature). 
[Q] proposed a model with periodically collapsing bubbles in which the bubble component of the 
price follows an exponential explosive path (the price being multiplied by at = a > 1) with 
probability it and collapses to zero (the price being multiplied by a t = 0) with probability 1 — tt. It 
is clear that, in this model, a bubble has an exponential distribution of lifetimes with a finite average 
lifetime vr/(l — tt). Bubbles are thus transient phenomena. In order to allow for the start of new 
bubble after the collapse, a stochastic zero mean normally distributed component b t is added to the 
systematic part of B t . 



Recently, pQ ] studied the implications of the bubble models for the unconditional distribution 
of prices, price changes and returns resulting from a more general discrete-time formulation of 
rational speculative bubbles: 

B t+ x = a t B t + b t , (1) 

where at can take arbitrary values and are i.i.d. random variables drawn from some non-degenerate 
probability density function (pdf) P a (a). The model can also be generalized by considering non- 
normal realizations of bt with distribution Pb(b) with E[b t ] = 0. In addition, the additive term 



is one way among many (see Q26| ] for alternatives) to nucleate the bubble from a non-zero value. 
Indeed, most of the time, a t < 1 and B t fluctuates with a scale set by bt- When an amplification 
phase starts up, B t increases exponentially from this initial value and the additive contributions of 
succeeding bt becomes subdominant in the exponential growth. In ([[]), B t denotes the difference 
between the observed price and the fundamental price which defines the bubble component. The 
term "bubble" refers to the regimes when B t explodes exponentially under the action of successive 
multiplications by factor at, at+i , ... with a majority of them larger than 1 but different, thus adding 
an additional stochastic component to the standard model of [Qj. 

Denoting E[.] the expectation operator, provided E[ln a] < (stationarity condition) and if there 
is a number fi such that < E[|b| M ] < +oo, such that 

E[|a|"] = 1 (2) 

and such that E[|a| M In |a|] < +oo, then the tail of the distribution of B is asymptotically (for large 
£>'s) a power law [ |l6| , ^] 

P B (B)dBn-^dB , (3) 

with an exponent fi given by the real positive solution of (||). Rational expectations require in 
addition that the bubble component in asset prices obeys 

<5E[B t+ i] = B t (4) 

where 5 is the discount factor < 1. This implies 

E[a] = 1/5 > 1 . (5) 



Since the function E[|a| M ] is upward convex, p0| ] showed that this automatically enforce fj, < 1. It is 
easy to show that the distribution of price differences has the same power law tail with the exponent 



fx < 1 and the distribution of returns is dominated by the same power-law over an extended range 
of large returns [p0[]. Although power-law tails are a pervasive feature of empirical data, these 
characterizations are in strong disagreement with the usual empirical estimates which find [i in the 
range 3 — 5 [34, 18, 24, 12, |i0[ ]. p0| ] concluded that exogenous rational bubbles are thus hardly 



reconcilable with some of the stylized facts of financial data at a very elementary level. This result 
has been extended to multi-dimensional bubbles involving an arbitrary number of coupled assets 



[21] 



Here, we provide a simple and natural extension of model (|J) that removes this constraint fi < 1 
and thus the discrepancy with empirical analysis. The key to understanding intuitively how the 
result of [ |20| ] derives is to realize that the rational expectations condition (Q) means that the bubble 
price grows locally exponentially with a instantaneous growth rate equal to the discount rate 

r s = - In 5 > . (6) 

Since the stochastic auto-regressive equation ([]]) describes a stationary process when the condition 
E[lna] < is satisfied, the only possibility to reconcile the non-stationary exponential growth with 
the stationary regression is that the distribution has no mean, which indeed occurs only for \i < 1. 
In practice for any finite time intervals, the absence of a mean implies that averages grow with the 
size of the time window. This remark points to the remedy that we now present. 



2 Generalization and breakdown of the constraint 
2.1 Rational expectation bubble model 

For the following discussion, it is useful to recall that pricing of an asset under rational expectations 
theory is based on the two following hypothesis: the rationnality of the agents and the "no-free 
lunch" condition. 

Under the rationnal expectation condition, the best estimation of the price pt+i of an asset at 
time t + 1 viewed from time t is given by the expection of p t +i conditionned upon the knowledge 
of the filtration {Tt} (i-e. sum of all available information accumulated) up to time t : E[p t +i\ Tt\. 

The "no-free lunch" condition imposes that the expected returns of every assets are all equal un- 
der a given probability measure Q equivalent to the historical probability measure P. In particular, 
the expected return of each asset is equal to the return rg of the risk-free asset (which is assumed to 
exist), and thus the probability measure Q is named the risk neutral probability measure. 

Puting together these two conditions, one is led to the following valuation formula for the price 

Pt ■ 

Pt = 8 ■ EQ[p t+ i|F t ] + d t V{p t }t>o , (7) 

where d t is an exogeneous "dividend", and 5 = (1 + rs)" 1 is the discount factor. The first term in 
the r.h.s. quantifies the usual fact that something tomorrow is less valuable than today by a factor 
called the discount factor. Intuitively, the second term, the dividend, is added to express the fact 
that the expected price tomorrow has to be decreased by the dividend since the value before giving 
the dividend incorporates it in the pricing. 

The "forward" solution of ([7]) is well-known to be the fundamental price 

+oo 

p? = X;* i -EQ[d t+i |F t ]. (8) 

i=0 



It is straightforward to check by replacement that the sum of the forward solution ^ and of an 
arbitrary component B t 

Pt =p{ + B t , (9) 
where B t has to obey the single condition of being an arbitrary martingale: 

B t = tf-E Q |B t+ i|F t ] (10) 

is also the a solution of (f7|). In fact, it can be shown (jllj] that @ is the general solution of (f7|). 

Here, it is important to note that, in the framework of the Blanchard and Watson model, the 
speculative bubbles appear as a natural consequence of the valuation formula ([7|), i.e., of the no 
free-lunch condition and the rationnality of the agents. Thus, the concept of bubbles is not an 
addition to the theory, as sometimes believed, but is entirely embedded in it. 



2.2 Model and predictions 

Consider again B t as the price difference between observed price and the fundamental price. We 
propose the following extension of the model initially defined by ([[]) which consists in introducing 
an exponentially growing additive term b t — ► e rt b t with r > such that the dynamics of B t is 

B t+1 = a t B t + e rt b t , . (11) 

The justification for modifying the additive term into e rt bt is the following. As we have seen, the 
additive term represents the background of "normal" fluctuations around the fundamental price. As 
we recalled above in the definition of the standard RE model, the additive term ensures that bubbles 
can nucleate spontaneously. If the economy and the market exhibit a long-term growth rate r, their 
"normal" fluctuations have also to grow with the same growth rate in order to remain stationary 
in relative value. The factor e rt in the additive term e rt bt thus reflect nothing but the average 
exponential growth of the underlying economy. In this framework, for the model to make sense, 
we need to impose an additional condition E[ln a] < r that the average growth rate E[ln a] of B t is 
smaller than the growth rate r of the background economy. This condition ensures that most of time 
B t is a "normal" fluctuation with stationary average relative deviations from the fundamental price. 
Only intermittent exponential amplifications will create transient bubbles that will be significantly 
above this (exponentially growing) background. 

We impose again E[b t ] = and bt is a white noise process. The rational expectations (Q) leads, 



as before [20], to 



E[a] = (T 1 . (12) 
We introduce the reduced price variable A t such that 

B t = e rt A t . (13) 

Equation (|ll|) then reads 

A t+ i = a t e~ r At + e' r b t , (14) 

which is of the standard form ([]]) with stationary coefficients and obeys the usual conditions. Note 
that the rational expectations condition ( fOj ) translates into E[a t e~ r ] = 5~ 1 e~ r which can now be 
smaller than 1 (see below). 



A 



Note that the condition E[lna] < r ensures that E[ln(ae r )] < which is now the stationarity 



condition for the process A t defined by (14). The conditions < E[|e r b t | M ] < +oo (which is the 



same condition < E[|b t |^] < +oo as before) and the solution of 

E[|ae"T] = 1 (15) 

together with the constraint E[|ae _r | M In |ae~ r |] < +oo (which the same as E[|a| M In |a|] < +oo) 
leads to an asymptotic power law distribution for the reduced price variable A t of the form Pa(A) ~ 
C J 4/|^4| 1+At , where \i is the real positive solution of (|T5|). Note that the condition E[ln(ae _r )] < 
which is E[ln(a)] < r now allows for positive average growth rate of the product ata t _ia t _2---«20iOo- 

[ |2"(i| ] have demonstrated rigorously that, if A t is distributed with a power law distribution with 
exponent fi, then its increments A t — A t -\ are also distributed according to a power law with the 
same exponent \i. Here, we are interested in the distributions of the returns of the price which we 
now discuss. 



2.3 Distribution of returns 

In order for the model to be meaningful with respect to returns, we need to recognize that the 
observable market price is the sum of the bubble component B t and of a "fundamental" price p{ 

Pt =p{ + B t . (16) 

Then, 

Pt+l = Pt+i + atB t + e r % = a t p t + (e r - a t )p{ + e r % , (17) 

where we have assumed that the economy and the market exhibit a long-term growth rate r, i.e., 
p t+ i = e r p t , leading to 

p{= Po e rt . (18) 



Expression (17) allows us to make concrete the meaning of the additive term e b t previously 



introduced in (|11|). Indeed, we pointed out that the additive term represents the background of 



"normal" fluctuations around the fundamental price (|18[). Their "normal" fluctuations have also to 
grow with the same growth rate in order to remain stationary in relative value. This is shown by 
replacing p{ in (|l6|) by poe rt given in (|l8|), which leads to 



Pt = e rt (p + A t ) . (19) 



In addition, replacing p{ in (17) again by poe rt leads to 

pt+i = a tPt + e rt \p (e r - a t ) + b t ] . (20) 



The expression ( |20| ) has the same form as (pT|) with a different additive term [po(e r — at) + bt] 
replacing bt- The structure of this new additive term makes clear the origin of the factor e rt intro- 



duced in (11): as we said, it reflects nothing but the average exponential growth of the underlying 
economy. The contributions e rt b t are then nothing but the fluctuations around this average growth. 

This novel additive term brings in a new feature which turns out to be essential to describe 
realistic return time series: it has a non-fluctuating component which ensures that p t grows on 
the average exponentially with time. For positive bubbles A t > or for large fundamental price 
Po > \Af \, a return ln(pt+i/pt) can not be made artificially large by the approach of p t arbitrarily 
close to 0. In constrast, the pure bubble component B t can cross an arbitrary number of times, as 



shown in figure 1. Close to these crossings, the returns [Bt+i — B t ]/B t of the bubble component 
are very large since the denominator is close to 0. Since such approach of B t to occurs with a 
uniform probability in the vicinity of 0, the distribution of 1/B t is determined by the rule of change 
of probability density under a change of variable and is found to be a power law with exponent 
equal to 1: this is one of the possible mechanisms known to generate power law distributions 
known as the "power law change of variable close to the origin" [ |30| ] (see Chap. 14). In the context 
of price returns, this is an artifact as the observable price contains the additional contribution of the 
fundamental price and therefore cannot go close to zero. 
The observable return is 

D Pt+i-Pt Pt+i ~ Pt + B t+i - B t (Pt+i-Pt , B t+1 -B t \ ( A t+ i-A 
R t = = f = Xt — 1 —f h - f =Xt[r + 



where 



f _ At f l f I At I l 

pt p{ + B t V pI Pt J \ p ° 

(21) 

Pt 1 

X* = f = i , (a I ^ • W 

p{ + B t 1 + (At/Po) 

In order to derive the last equality in the right-hand-side of <^1\), we have used the definition of the 
return of the fundamental price (neglecting the small second order difference between e r — 1 and r) 
and we have combined (|l3|) and (18). Expression (^l|) shows that the distribution of returns R t of the 
observable prices is the same as that of the product of the random variable xt by r+ (At+i — A t ) /po. 
Now, the tail of the distribution of r + (A t+ i — A t ) /po is the same as the tail of the distribution of 
At+i — A t , which is a power law with exponent \i solution of (|l5|), as shown rigorously by [p0|]. 

It remains to show that the product of this variable r + (At+i — A t )/po by x has the same tail 
behavior as r + (At+i — A t ) /po itself. If r + (At+i — A t ) /po and x were independent, this would 
follow from results in [||] who demonstrates that for two independent random variables 4> an d X 
with Proba(|^>| > x) m cx~ K and E[x K+<E ] < oo for some e > 0, the random product <p% obeys 
Proba(|0x| > x) w E[x K ]x~ K . 

r + (At+i — A t )/po and x w& not independent as both contain a contribution from the same 
term A t . However, when A t « po, x is close to 1 and the previous result should hold. The impact 
of At in x becomes important when A t becomes comparable to p$. 



It then convenient to rewrite (21 ) using (22) as 



= r A t+1 - At = r {a t e~ r - l)A t + e~ r b t 

1 l + (V?o) Pk + M 1 + Wpo) Po + At ' { ' 

We can thus distinguish two regimes: 

• for not too large values of the reduced bubble term A t , specifically for A t < po, the denomi- 
nator po + A t changes more slowly than the numerator of the second term, so that the distri- 
bution of returns will be dominated by the variations of this numerator (a>te~ r — l)A t + e~ r b t 
and, hence, will follow approximately the same power-law as for A t , according to the results 
of [§. 

• For large bubbles, A t of the order of or greater than po, the situation changes, however: from 



(|23[), we see that when the reduced bubble term A t increases without bound, the first term 
r/(l + (At/po)) goes to while the second term becomes asymptotically ate~ r — 1. This 
leads to the existence of an absolute upper bound for the absolute value of the returns given 

by 

max|i?f| = max{|min(ate r — 1)|, |max(ate r — 1)|} . (24) 



To summarize, we expect that the distribution of returns will therefore follow a power-law with the 



same exponent \i as for A t , but with a finite cut-off given in equation (24). 

Consider the illustrative case where the multiplicative factors at are distributed according to a 
log-normal distribution such that E[lna] = lnao (where ao is thus the most probable value taken 
by at) and of variance a 2 . Then, 



E[|ae~T] =exp 



2^ 

-r/x + /i In a + \x — 



(25) 



Equating (|5|) to 1 to get \i according to equation (|i"5|) gives 



r-lna r-lna r - r s 

V = 2 5 = : = H : . (26) 

a z rs — in ao rg — in ao 

We have used the notation @ for the discount rate defined in terms of the discount factor. The 
second equality in ( p6| ) results from (|I2|) using E[a] = aoe' 7 I 2 . 

First, we retrieve the result [ 20 ] that a < 1 for the initial RE model (|l|) for which r = and 



In ao < 0. However, as soon as r > rs = — In 5, we get 

H > 1 , (27) 



and fx can take arbitrary values. Technically, this results fundamentally from the structure of ( 1 1 ) in 



which the additive noise grows exponentially to mimick the growth of the bubble which alleviates 



the bound \i < 1. Note that r does not need to be large for the result (27) to hold. Take for instance 
an annualized discount rate r v s = 2%, an annualized return r y = 4% and ao = 1.01. Expression 
predicts /x = 3, which is compatible with empirical data. 



2.4 Implications of the condition r > r s 

When the price fluctuations associated with bubbles grow on average with the mean market return 
r, we find that the exponent of the power law tail of the returns is no more bounded by 1 as soon as 
r is larger than the discount rate r<j and can take essentially arbitrary values. As can be seen from 
equation (|8|), this condition r > rs corresponds to the unsolved regime in fundamental valuation 



theory where the forward valuation solution (||) loses its meaning, as discussed recently in p8|]. 
Indeed, changing d t +i into dt+ie ri to account for the growth of dividend associated with the growth 
of the fundamental price, the sum in (||) then behaves as J2t^o [^ e T which diverges for 5e r > 1, 
i.e., r > rs (neglecting the difference between e r and 1 + r). 

Two attitudes can be taken with respect to the condition r > rs, entailing two different modeling 
strategies. 

2.4.1 The "hard-line" attitude 

The "hard-line' attitude (summarized from private exchanges with T. Lux) can be described as 
follows. The assumption r > r$ corresponds to abandon the background of the rational bubble 
theory which is the rational valuation formula which fails. Once this starting point is abandoned, 
there is no theory of rational bubbles anymore as the transition from (|h to @ with (|I0|) is no 



more defined. The question of whether the results of Q2O0 are saved or changed then becomes 
irrelevant. In other words, giving up the fundamental valuation formula, the argument goes, would 
mean that we give up "rationality" (in the very limited and special sense of the rational valuation 



formula and similar theories in economics). We would then find ourselves within the realm of not- 
fully-rational behavior which implies that we have already given up also the possibility of having 
"rational speculative bubbles". This does not preclude any kind of near-rational bubbles like those 
in the agents-based models (see for instance [19, 23, 13, 171 31]) but the strictly rational bubbles 



cease to exist. To sum up the hard-line attitude, the theory of RE bubbles is of no interest anymore 
when one gives up the hypothesis of rational valuation because there is then readily available a 
universe of possible alternative bubble theories. 

2.4.2 The "intrinsic RE bubble" model 

The other attitude advocated here and elsewhere 



[ |14] , \l3j , p2| , [33]] decouples the rational valuation formula from the RE bubble model, which we term 
the "intrinsic RE bubble" model. It is true that, in the context of rational pricing based on the flow 
of dividend with infinite time horizon, an exponential growth of the price is associated with the 
same exponential growth of the dividend. This leads to an inconsistency as the valuation formula 
leads to a diverging price if the growth rate is larger than the discount rate. As discussed in [E^J, the 
price is of course not infinite, only the rational valuation formula is of no use to determine it. The 
price has to be defined by processes other than the dividend valuation formula. 

• From analogies with statistical physics and with the theory of bifurcations and their normal 



forms, in which similar situations occurs, I proposed in [28] a scenario according to which 
this regime may be associated with a spontaneous symmetry breaking phase corresponding 
to a spontaneous valuation in absence of dividends. 

Another approach is to realize that the fundamental valuation formula is a statement of equi- 
librium. It is thus natural to interpret the regime r > r$ where the fundamental valuation 
formula fails as intrinsically dynamical: the fundamental price becomes determined by dy- 
namical processes other than just the flow of dividends. 

It is possible to develop models of general equilibrium [Q] in which one can treat the growth 
of the economy reflected in the stock market as the accumulation of composite capital, plus 
a change of available technology that affects both the old and new capital. Capital's market 
value grows as people add new units of capital, or as existing units grow in value as they 
become more effective. If K t is the composite capital or market value at time t, r t is the 
return and q is consumption, then K t+ i = r t K t — ct- In the long run, consumption seems 
to grow at the same rate as wealth, income and other measures of available resources. Thus 
q is roughly proportional to K t , i.e., a = 7t-Kt with j t < Tt as people in general con- 
sume less that the expected return on composite capital, leading to expected growth. In this 
framework [0], dividends are seen as depending on current and past prices rather than prices 
as depending on expected future dividends as in the standard valuation formula. It thus is 
perfectly possible to model a market growing at an average rate r equal to the average of 
Tt — It which can be larger than the risk-free rate or discounting rate or inflation rate. This 
does not invalidate the assumption of rational expectations, which simply states that people 
make economic decisions in a way that tends to take into account all available information 
bearing significantly on the future consequences of their decisions. We can still form rational 
expectations from an intelligent appraisal of circumstances, though the process behind such 
circumstances may be hard to discern. An economic system doing the most efficient possi- 
ble job of reading the information being reflected in price signals will still experience some 



irreducible business cycle swings and bubbles. This results from the fact that the economic 
process contains inherent mechanisms that convert random shocks on prices into a more per- 
sistent, short-term misreading of changing profit opportunities, in other words the rational 
interpretation of noise can stimulate a cumulative swing in output that will continue until the 
misreading is realized and retrenchment sets in. Random shocks to prices and markets are 
always with us and rational bubbles are thus naturally occurring market phases. Burke [||] has 
also established, by dropping continuity and weakening utility representation, that commod- 
ity prices and consumptions can approach approximate equilibrium to within any practical 
tolerance after dropping the standard general-equilibrium assumption that preference orders 
discount future consumption faster than the economy grows. 

• Another fix is based on the understanding that the divergence of the sum in (|8]) stems from the 
fact that the investor assumes an infinite time horizon; in practice, we would discount future 
cash flows only up to a maximum time, say ten or thirty years at most. As a consequence, 
the rational valuation formula has to be truncated to a finite number of terms and provides a 
finite price. 

Therefore, the apparent contradiction is resolved by abandoning the pricing based on dividends 
discounted over an infinite time horizon and replacing it by an exogeneous dynamical process. In 
other words, we assume here that the fundamental price is growing exponentially, without reference 
to dividends. We thus propose not to "throw away the baby with the bath", i.e., to decouple the 
fundamental pricing mechanism from the RE bubbles. Indeed, investors use various methods for 
estimating the fundamental price and the rational valuation formula is only one among an arsenal 
of techniques that, in real life situations, are usually combined. 

The point of view advocated here is that it is reasonable to keep the RE bubble model (JTo|), 
independently of what determines the specific underlying fundamental price. This is actually the 
point of view initially introduced by ^ and further developed in [14, |J, 33], where the bubble 



price may follow essentially arbitrary trajectories as long as there is a jump process (crash) with a 
hazard rate ensuring the no-arbitrage condition. 

In a way, the choice between the "hard-line" attitude and the "intrinsic RE bubble" model 
is almost a matter of philosophy with respect to the meaning of models. Ultimately, a model is 
useful if it helps thinking about the problem and if it allows one to formulate new questions. It 
seems interesting to use the modelling strategy in which the RE condition is relaxed partially by 
abandoning only the rational valuation formula. Following this view point, we have the following 
logical sequence: 

1 . RE bubbles give an exponent \i > 1 for the distribution of returns only when the fundamental 
price is growing faster than the discount rate; 

2. when the fundamental price is growing faster than the discount rate, the rational valuation 
formula breaks down; 

3. empirical studies show that fj, > 1 in real financial time series; 

4. as a consequence, since the distribution of returns are usually constructed over rather long 
time histories, it is sufficient that the fundamental valuation departs from the dividend formula 
at some times to explain item 3. 

This provides a parsimonious model for the empirical distribution of returns. 



2.5 Numerical simulations 



We now present numerical simulations that illustrate these results. We do not pretend to capture 
reality acurately but show nevertheless that the stylized facts of empirical data are recovered by 
this simple model. In order to make a precise comparison with empirical data, one should specify 
the possible shapes of the distributions P a (fl) and P&(6) and one should also add the fundamental 
price to the bubble price B t since only the sum is observable. Tests of the model would thus be a 
joigned test of the relevance of the RE model together with an assumption on the dynamics of the 
fundamental price. 

Figure [j] shows a typical trajectory of the price B t of a bubble generated by equation (11) with 



r = 0.33%. This choice is such that one time step can be approximately interpreted as one month, 
since the compounded return over 12 time steps gives a realistic yearly return of 4% (1.0033 12 = 
1.04), as can be seen in figure 2. 

The multiplicative factors at are generated from the formula 



at = ao ex.p[ar] t ] , where ao = 1.001 and a = 0.0374 , 



(28) 



and the rjt's are centered Gaussian random variables with unit variance. With this parameterization, 
at are log-normal random variables of variance a 2 with E[lna] = lnao = 0.001 and E[a] = 
a e CT2 / 2 = 1/5 = 1.0017. Therefore, 6 = 0.9983 and r s defined by @ is r$ = 0.17%. With the 
values r = 0.33%, rs = 0.17% and lnao = 0.1%, expression ( [26] ) predicts an exponent fi = 3.3. 

The additive term b t is taken uniformely distributed in the interval [—0.05; +0.05]. Note that b t 
sets the scale of B{t). The bubble price B t stays on average at 0, as seen from (JTTJ) which predicts 
E[B t ] = 0. However, like a correlated random walk with fat tails, it wanders around, with an 
appearance quite reminiscent of active and volatile markets. 

As discussed in the previous section, in order to make a meaningful comparison with empirical 
data, it is necessary to add a fundamental component to the bubble to form an observable price. 



To remain as simple as possible, we follow the specification (16) with (18), with po = 1- Figure 
H shows the time dynamics of B t + e rt , where B t is the same as shown in figure [I] and we have 
assumed that the fundamental price follows a deterministic growth at the "annualized" rate of 4% 
(corresponding to the monthly rate r = 0.33%) according to (|TJ). The trajectory of the Dow 
Jones Industrial Average (DJIA) extrapolated back from 1790 till sept. 2000 is also shown as 
the thin line. The Dow Jones index was constructed by The Foundation for the Study of Cycles 



( http : j /www. cycles. org /cycles. htm ). It is striking to observe how the simple RE bubble model 
added to a simple exponentially growing fundamental price can capture the large scale variability 
of the DJIA. 

Figure |3| shows a double logarithmic scale representation of the complementary cumulative 
distribution of the "monthly" returns R t defined in (^|), constructed from the synthetic time series 
shown in figure ^[ As predicted, the distribution is well-described by the asymptotic power law 



with an exponent in agreement with the prediction fi 3.3 given by the equations < \15\ ) and 
and shown as the straight line. Note that the expected cut-off at the maximum return is not observed 
as the returns have not yet explored the large values of the order of 1 in the finite-time series. 

For comparison, figure |] shows the distribution of positive and negative returns of the Dow 
Jones Industrial Average price over the 20th century. The tails are very well quantified by power 
laws with exponents respectively equal to fi + = 2.9 ± 0.3 and /u_ = 2.4 ± 0.3. The error bars 
are estimated from the theory of maximum likelihood applied to the Hill estimator which predict 
5/j,//j, = 1/ y/~N, where N is the number of data in the power law tail. Using N m 100, we find Sfi « 
0.3. The data does not reject the hypothesis that fi + = /i_ 2.7±0.3. These values are compatible 



with those previously reported in the literature on smaller time scales [ p4[ , |18| , [24] , |12| , |i0( |. The 
comparison between figures [| and ^| suggests that the RE bubble model (11) goes a long way into 
capturing the structure of real market prices. 



3 Concluding remarks 

In the numerical example given in the previous section, the stationary process A t with zero mean 
and finite variance may in principle eventually reach the value — po with very low but not zero 
probability, at which the observable price vanishes. Such large negative bubbles will thus produce 
artifacts in the distribution of returns which, as we demonstrated for the returns of B t , are controlled 
by the excursions of p t close to 0. The general formulation of the RE bubbles allows in principle 
for both positive and negative bubbles with the sign of each new bubble depending on that of the 
additive stochastic term in its inception period. However, it is well known that there are conceptual 
problems with negative bubbles [§]. Therefore, we deliberately confined ourselves to bubbles of 
amplitudes smaller than the fundamental price. Numerical tests using only positive bubbles, ob- 
tained by imposing a strickly positive additive term b t , confirm the results of figures 2 and 3. We 
thus emphasize that none of our results hinges on this slight modification of the model and both 
modeling strategies are equally viable: consider only the bubbles that do not become too negative 
so that the total observable price does not become too small or construct only positive bubbles. This 
conclusion has also been checked extensively by ^Ofl. 

While investors enjoy getting a larger return r, this comes usually at the price of increasing risks. 
Here, the situation is different because a higher market return leads to a thinner tail of the return 



distribution since the exponent [i, for instance given by fl26[), increases with r. Hence, increasing 
the average market return r decreases the extreme risks. In this context, it is interesting to bring 
into focus the long-standing paradox that the Dow Jones average index has been argued to exhibit 



an anomalously large return, averaging 6% per year over the 1889-1978 period [g2J, which cannot 
be explained by any reasonable risk aversion coefficient. What our study shows is that the excess 
remuneration r — r$ has an unexpected "good" consequence in decreasing drastically the large risks 
of the market by increasing the exponent of the asymptotic power law distribution. Following recent 



ideas in the theory of complex systems (see for instance [ pQ| ] and references therein), we argue 
that the market has self-organized such that the excess remuneration has reached a level which 
compensate for the huge risks associated with the intermittent bubbles created by investors. We 



conjecture that the presently observed value of the exponent \i in the range 3 — 5 [ p4[ , |18| , [24], [12], [10|] 
and the anomalous returns of the market over long period of times are the tools that the system has 
developed to tame the fat tails that bubbles tend to create. Extending the present model to derive a 
dynamics on r would allow us to understand how this organization may proceed. 

Acknowledgements: I am grateful to T. Lux and Y. Malevergne for stimulating discussions. 
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Figure 1: Bubble price B t generated by equation (Mil) with r = 0.33%. The multiplicative factors 



at are generated from the formula (28). a t are thus log-normal random variables with E[lna] = 
lnao = 0.1% and E[a] = 1.0017. The additive term bt is uniformely distributed in the interval 
[—0.05; +0.05]. With these parameters, one time step can be interpreted as approximately one 
month. The time interval shown covers thus approximately 80 model-years. 
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Figure 2: The thick line shows the time dynamics of B t + p t (bubble + fundamental price), 
where we have assumed that the fundamental price p{ = e rt follows a deterministic growth at the 
"annualized" rate of 4% corresponding to the value r = 0.33%. The trajectory of the Dow Jones 
Industrial Average (DJIA) extrapolated back from 1790 till 2000 is also shown as the thin line. 
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Figure 3: Double logarithmic scale representation of the complementary cumulative distribution 
of the "monthly" returns R t defined in ( |2l| ) of the synthetic total price shown in figures ||. The 
continuous (resp. dashed) line corresponds to the positive (resp. negative) returns. The distribution 
is well-described by an asymptotic power law with an exponent in agreement with the prediction 
fi « 3.3 given by the equations (JT3|) and (E3) and shown as the straight line. The small differences 



between the predicted slope and the numerically generated ones are within the error bar of ±0.3 
obtained from a standard maximum likelihood Hill estimation. 
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Figure 4: Double logarithmic scale representation of the complementary cumulative distribution 
of the monthly returns hipt+i/pt of Dow Jones Industrial Average price p t from Jan. 1900 till 
Sept. 2000. The +'s (resp empty circles) correspond to positive (resp. negative) returns. The 
two dashes straight lines give the best fits to a power law for the largest 100 positive and negative 
monthly returns. For positive (resp. negative) returns, we find an exponent /i + = 2.9 ± 0.3 (resp. 
/i_ = 2.4 ± 0.3). This is in agreement with previous empirical works [ |34] , 18, 24, [12|, 10]. 



1 -7 



